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1.  Introduction 

The  purpose  of  this  paper  is  to  put  together  in  a  concrete 
form  the  mathematical  theory  of  an  inverse  scattering  problem  which 
was  originally  formulated  in  quantum  mechanical  terms  but  which  is 
related  to  the  problem  of  determining  the  variation  of  electron 
density  in  an  ionized  gas  such  as  the  earth's  ionosphere  by  means  of 
radio  sounding  experiments.  The  theory  described  here  was  developed 
by  H.  E.  Moses  and  the  present  writer  and  was  originally  presented 
I  rather  generally  and  abstractly. In  the  present  article  the 

definitions  and  proofs  of  theorems  have  been  revised  in  an  attempt  to 
make  the  analysis  understandable  to  a  reader  who  may  not  be  familiar 
with  the  terminology  of  quantvim  mechanics  or  of  the  associated  Hilbert 
space  theory.  While  the  usual  terminology  employed  by  physicists  has 
not  been  discarded  completely  here,  definitions  and  relations  are  given 
in  concrete  terms  whenever  possible,  and  abstract  relations  are  used 
for  convenience  of  notation  rather  than  for  the  purpose  of  obtaining 
great  generality. 

The  physical  problem  at  which  the  present  work  is  aimed  is  to 
find  the  variation  of  electron  density  in  a  weakly  ionized  gas  from 
a  knowledge  of  the  scattering  amplitude  resulting  from  the  incidence  of 
a  plane  electromagnetic  wave.  It  is  assumed  that  the  relevant  properties 
of  the  medium  vary  slowly  enough  so  that  a  scalar  wave  describes  the 
scattering  phenomena  with  sufficient  accuracy.  It  is  also  assumed  that 
j  losses,  ordinarily  attributed  to  the  collisions  of  the  electrons  with 
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heavy  particles  may  be  neglected.  The  effect  of  an  external  magnetic 
field  which  would  cause  anisotropy  is  neglected  as  well,  so  that  the 
scattering  is  assumed  to  be  isotropic  as  well  as  lossless. 

A  steady-state  scalar  field  component  u(x,y,z)  under  these 
conditions  will  satisfy  a  differential  equation 


(A) 


An  +  Jk  -  V(x,y,z)l  u  =  0, 


having  the  same  form  as  the  time-independent  Schrodinger  equation  of 

[3] 

quantum  mechanics.  The  quantity  k  is  the  wave  number  and  the  function 
V'(x,y,z)  depends  upon  the  variation  of  electron  density  in  the  medium; 
in  fact,  V(x,y,z)  is  the  plasma  frequency  divided  by  the  square  of  the 
velocity  of  light  in  vacuum. 

The  scattering  problem  associated  with  such  a  differential 
equation  would  be  to  calculate  the  scattering  amplitude  T(x,y,z,k) 
defined  by  the  asymptotic  form  of  a  solution  u  for  which  the  incident 
wave  is  plane  and  -uniform  and  propagates  in  the  direction  of  the  vector 
k.  That  is,  if 


u.  =  e 
inc. 


ik'X 


where  x  is  the  radius  vector  whose  components  are  x,y,z,  and  at  large 
distances 


(B) 


ik'X  ,  _/■  ,  \  e 

u  ~  e - +  T(x,k)  — 


ikx 


where  k  and  x  are  the  magnitudes  of  the  vectors  k  and  x,  the  inverse 
scattering  problem  would  be  to  determine  the  function  V(x,y,z)  in  (A), 
given  the  scattering  amplitude  T(x,k),  or  some  part  of  it. 

At  the  present  time  the  inverse  scattering  problem  for  the 
system  described  precisely  by  equation  (A),  in  three  dimensions,  iS 
still  open.  In  this  paper  we  shall  consider,  instead  of  (A),  the 
more  general  integro-differential  equation 

(C)  Au  +  k^u  -  J v(x,x' )u(x' ,k)dx'  =  0, 

where  the  integration,  as  indicated  by  the  differential  dx'  is  over 
all  three-dimensional  x'  space.  ¥e  assume  that  the  integral  operator 
whose  kernel  is  v(x,x' )  is  such  that  the  asymptotic  relation  (B) 
still  holds  for  some  solution  of  (C).  Physically,  the  differential 
equation  (C)  would  govern  electrcxnagnetic  wave  propagation  in  a 
somewhat  artificial  plasma  of  a  more  general  type  than  that 
characterized  by  (A).  However,  we  allow  the  kernel  v(x,x'  )  to  be 
a  distribution,  and  in  partictdar  it  could  have  the  form 

(D)  v(x,x')  =  V(x)  &(x-x')> 

whereupon  the  integro-differential  equation  (C)  would  be  indentical 
with  the  differential  equation  (A). 

From  the  point  of  view  of  the  inverse  scattering  problem  we 
woiild  have  to  discover  a  method  of  solution  which  would  automatically 
impose  the  condition  (D)  if  we  wished  to  solve  the  original  problem 
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associated  with  (A).  Instead,  what  we  actually  do  here  Is  to  impose 
conditions  that  guarantee  a  unique  solution  of  the  inverse  problem 
associated  with  (C)  but  which  do  not  guarantee  the  form  (D)  in  general, 
although  the  possibility  is  not  excluded  that  for  some  given  scattering 
operators  the  solution  v(x,x' )  will  have  the  form  (D)  and  hence  will 
apply  to  the  problem  associated  with  (A). 

In  section  2  of  this  article  various  conventions  which  are  used 
throughout  are  listed  for  convenience  and  are  referred  to  specifically 
later  on  whenever  it  seems  appropriate  to  do  so. 

In  section  3  an  integral  transform  theorem  analogous  to  the 
Fourier  Integral  theorem  and  involving  arbitrary  solutions  u(x,k)  of 
(C)  is  presented.  This  theorem  requires  the  use  of  a  certain  weight 
function  w(k,k' )  needed  to  normalize  the  particxilar  function  ^(x^k) 
"which  appears  in  each  Instance  of  the  theorera»  '  We  calculate  'w(k,k' ) 
for  the  case  in  which  u(2C,k)  is  the  ^function  satisfying  the  asymptotic 
condition  (B).  The  discussion  should  be  of  seme  interest  by  itself 
for  other  more  usvial  applications  such  as  the  formation  of  wave  packets 
in  the  medium  determined  by  (C)  (or,  in  particular,  by  (A))  associated 
with  normalized  scattered  plane  waves. 

Section  (4)  deals  with  an  inverse  scattering  problem  for  (C) 
wherein  conditions  sufficient  to  guarantee  uniqueness  of  the  solution 
are  imposed.  It  is  hoped  that  the  same  or  a  similar  formal  setup  can 
be  used  in  the  future  to  attack  the  inverse  scattering  problem  asso¬ 
ciated  directly  with  (A). 
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2.  Notation 

In  this  section  we  list  a  set  of  conventiais  for  easy  reference 
later  on.  The  notation  described  here  will  be  used  consistently 
throughout  the  rest  of  this  article. 

a)  An  underlined  letter  such  as  x  will  represent  a  three-dimensional 
vector  (x^^XgjX^).  The  same  letter  x,  hot  underlined,  will  represent 


/  2  2  2- 

the  magnitude  v  x^  +  x^+x^  of  that  vector. 


3 

b)  Unless  there  are  explicit  limits,  an  integral  will  be  over  the 

entire  dcxiiain  of  the  variable  occurring  in  the  differential.  Thus.' 

00  00  00  00  .  ' 

J' f(x)dx  means  J' f(x)(ix:,  while  J' f(x)dx  means 


f ( x^ , Xg , )dx^dx2dx2 . 


•00  -00  -00 


Mxiltiple  integrals  will  be  designated  by  the  nature  of  the  differentials 
which  occur  in  them,  and  a  single  integral  sign  will  always  be  used. 

For  example,  J'  is  sl  six-dimensional  integral  oyer  the  entire 

domain  of  the  vectors  x  and  y. 

c)  We  shall  be  concerned  with  linear  operators  in  the  form  of 
integral  operators  over  a  suitable  space  of  functions.  The  kernel  of 
such  an  operator  will  be  denoted  by  a  lower  case  letter.  When  we 
wish  to  speak  of  the  operator  in  the  abstract,  the  corresponding 
capital  letter  will  be  used.  Thus,  Uf  means  the  abstract  operator  U 
^acting  on  the  abstract  linear  element  f,  and  the  concrete  representer 
of  U  will  be  J' u(x,x'  )f(x'  )<ix'* 
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d)  The  kernels  allowed  for  integral  operators  can  be  distri¬ 
butions  as  well  as  smooth  functions.  Thus,  in  particifLar,  the 
identity  operator  I  is  represented  by  a  kernel  which  is  a  delta  function 


If  =  J S(x  -x'  )f(x'  )dx'  =  f(x)  =  f. 


The  Laplaclan  operator  A  can  be  represented  by  an  integral  operator 
whose  kernel  is  given  by  A6(x- x'  )  =  A*6(x-  x'  )>  wherein  the  primed 
differential  operator  is  taken  with  respect  to  the  primed  variables, 
and  the  imprlmed  operator  with  respect  to  the  unprlmed  variables.  .For  , 
example,  we  have 


Af  s  Af  (2) 


I 


f(2'  )dj£'. 


e)  We  shall  use  the  notation  U*’  for  the  Hermit ian  adjoint  of  U, 
but  u*(x^x'  )  will  mean  the  complex  conjugate  of  the  function  u(x,x'  )., 
Hence  the  following  are  equivalent: 


U*f  =  J' u*(x' ,x)f  (x' )dx'  . 


f )  Instead  of  functions  of  2  the  Fourier  transforms  of  the 
functions  will  sometimes  be  used  to  represent  the  same  linear  space. 
The  new  functions,  the  transforms,  will  depend  on  vectors  k  instead  of 
X  and  will  be  indicated  by  a  circumflex.  Thus 


fsf(x),  f  =  f(k) 


1  r  ik-x  V- 

^  /  e - f(x)dx. 
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In  the  case  of  an  integral  operator  kernel  we  shall  also  write 

'^'3/2  J '^(2>2£’)  dx'. 

Because  of  the  identity 

— .  f  e^-  ^  ^  dk  =  &(x  --x'  ) 

i2Tty,.^ 

we  have  .r 

.  ■  '  J' 

g)  The  Laplacian  operator.  A  which  acts  on  the  functions  off  (x) 
is  represented  in  k-space,  the  space  of  the  Fourier  transforms.  f(k),  '/ 

by  the  multiplicative  operator  -k  (cf.  a),  that  is,  by  the  integral 
operator  whose  kernel  is  rk  5(k-k'  ).  .  ..Thus,. 

Af(x)  =  -  ■  f(k)dk  =.-  --  ■■■  7 ■■  Tkhik-  k'  )e"^- -  f(k' )dk'dk. 

.Hence  we  can  write  , 

;  AH-k^=  -k^6(k-kV) 

in  k-space.  Because  of  the  relation 

f(x)S(x)  =  f(0)S(x) 
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any  operator  whose  kernel  has  the  form 


r  =  r(k  -k' )  &(k^-  k'^) 


commutes  with  the  Laplacian  A.  .  That  is. 


^  =  FA, 


for  we  have 


;  .  .  -k^r(k  -k'  * >k'^r(k-,k  k'^). 


3.'  Eigenfunctions  ' 

Solutions  u  (x,k)  of  theXdifferential  equation 

O  ' 

(1),  -A  u^(x,k)  “  k^u^(x,k)..  . 


will  be  called  eigenfunctions  of  the  operator  -A.  Each  solution  , 

•  u^(x,k)  will  be  regarded  as-  the  kernel  of  an  integral  operator  which 
Is  a  representer  of  some  abstract  operator  tj^.  This  integral  operator 

acts  on  the  space  of  functions  f(k)  which  are  transfoiroed  by  it  into 

functions  f(x)  belonging  to  the  x-representatlon.  Thus,  we  shall  have 
relations  of  the  type 

A  special  case  of  this  is  the  Fourier  transform  for  which 


/  .  ^  1  -ik*x 

u  (2>k)  =  o/p  e  -- 


(cf.  2.f). 
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Similarly,  we  shall  be  concerned  •vrtth  operators  U  and 
corresponding  eigenfunctions  u(x,k)  which  are  solutions  of  the 
integro-differential  equation 


(2)  7  ■v:(x,,x' )u(x' ,k)dx'  =  k  u(x,k). 


In  operator  form  (2.)  is 


-AU  .+  ':VU  =  k  U. 


Here  again  u(x,k)  is  regarded  as;  the  - kernel  of  an  integral  operator 
which  tr'ansfoimis  ■  a  function  of  k  into  a,  function  of  x*  In  (2)  the 
operator  V  whose  kernel  is  v(x,x'  )  may!  in  Special  cases  be  a  multi¬ 
plicative  operator;  i.e. , .  v(x,x' )  may-have  the  form.  Cp(x)6(x-x*  )  .- 
.■■'(cf.  2.'d).  .  Then  (2)  will,  have  , the -usual  form'  of  the  Schrodinger '  eq-iiatiop 

■■'■y/'V  It.  is  well,  known' .that  the  operator— A  ..is 'self-adjoint  and.  has 
■a ,,  continuous  spectrum  consisting,;'Of  'zero' and  all  positive  real  hmber's'.' 
'when,  acting  on  a  space,,  e.g. ,  of  square  integrable  functions  f(x).  'We  . 
shall  ass-um'e  that  this  is  true  of  .'-A  and  in  addition  of  -A  +  'V  without 
worrying  about  particular  conditions  :on  the  "perturbation"  Y  necessary 
to  guarantee  this  property  or,,  for,  that  matter,  about  the  precise 
nature,  of  the  space  of  functions  f(x).  t! 


.  .  We  have  given  ..an',  operator-  interpretation  of  the  eigenfunctions 
of  (1)  and  (2)  as .transformations  fran  k  space  to  x  space.  The  more 
concrete  interpretation  is  also  useftil,  especially  in  applications,  in 
which  we  regard  them  as  providing  integral  transform  theorems  analogous 
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to  the  Fourier  Integral  theorem.  We  can  then  use  them,  for  example, 
to  describe  the  transient  response*  to  the  system  described  by  (l)  or 
(2).  Then,  the  Fourier  Integral  theorem  provides 


(3)  ./  .f(x)  =  ~“"T7p  r  ^ 

■  ..  ..  .  J,  -  \, 


tiK)dk, 


f(k)  =  — 


If  .we  set 


then. Q)  and  (4)  Imply  the  operator  relation  , 

where  i  is  the  identity.  ■  ,  The  relations  (5)  mean  that  is  unitary^ 

'  Similarly  If- we  define  the  particular  eigenfuiiction  u^(x,k) 

Of  :(2)  as  the  unique^  solution  of  the;  inte^al  equation  , 

then  we  can  prove  an  analogous  transform  theorem  in  terms  of  u_^(x,k) 

The  transient  response  would  be  given  by  multiplying  a  function  f(k) 
by  e^'*’^u(x,k),  where  w  =ck,  and  integrating  oyer  w. 

^The  solution  of  the  integral  equation  (6)  is  mique  because  of  our 
assumptions  about  the  operator  -A  +  V. 
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and  u*(x,k).  We  can  prove: 


if 

(7) 

then 


f(x)  =  J  ^i+(x,k)(p(k)dJk, 


(8)  ;  ■■  ^•V'>(k) 

Relations  (7)  arid;  (8)  imply  . the  i-el^^  ... 

:'(9)  'u^*  =\^ 

for  .the  operator  U^. corresponding  to  u^(x,k)i  Thus,  is  unitary. ■ 

Similar  results  hold' for  the  eigenfunction  u_(2,k)  and  the 
corresponding  operator  U_,  where  ■•u_(x,k)  is  defined  as  the  solution  of 
the .  integral  equation 


ikix-x' 


T  •  •  4 1^'  V  •  T  “*liv  I  X*X  I 

-  ki  J  v(a',x")u.(s%s)te'4" 

y.We  shall  prove  the  assertion  we  have  made  concerning  the  eigen- 
factions  u+(x,k)  by  showing  that,  for  example, 

(11)  •  y'.u*(x,k  )u^(x»k’  )<3x  =  S(k  -  k' ). 


Since  we  have  assumed  that  the  spectra  of  -A  and  -A  +  V  are  the  same 
the  self-adjointness  of  these  operators  will  imply  that  has  a  right 
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inverse  as  well  as  a  left  one.  Thus,  because  of  (11)  we  shall  have 


* 


=  I- 


The  arguments  required  to  prove  a  similar  theorem  in  terms  of  U  are 
canpletely  analogous,  and  we  shall,  have  finally 


(12) 


uX  “  uX  =  1/ 


or 


^i+(x,k)u+(^,^')d2  =  S(li.-k'), 


(13) 


“  ^(2£“  2£'  )• 


i^oof  of  (ll)t  •  Before  beginning  the  argment  it  will  be  useful 
to  define.  ' 


■v  ■  y  . 


t(a,]i)  =  /  v(2,3s'  )u^(i'  ,^)da' . 


We  can  now.  vrlte  in  place  of  (6) 


.  •  1  1  r  ik(x-2£'  I 


(2k) 


We  shall  eJ.so  need  the  well-known  identity 


(15) 


^  3±ikU-a' 
kit  lx  -x' 


lim 
e  -»  0 


ip.  (jc-x' ) 


e. 

(2jt)'^  ^  k^-  p^  ±i€ 


dp. 
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which  is  proved  in  appendix  I.  The  relation  (].5)  has  the  meaning  of 
a  distribution;  i.e.,  the  limit  is  to  be  taken  after  an  integration 
over  the  variable  x '  : 


We  shall  also  require  the  identity 


.(i6) 


(a±ie)(b±ie)  “  (a+b)±2ie 


1 


a±i€  b±ie 


which  follows  frati  a  simple  straightforward  a,lgebraic  manipulation. 


In  addition  we  shall  need.-.  . 


(17) 


— ilrCx'  ,k''  }dx' 

ie 


and 


"1  ■■  r I  '■ ''i  ■  ,  Tw .r  e^”  ■  */■ "  ' 

-  S'  f  \  "  ^ 

..  y-...  ■■  .  X  .-x'  ■  .v^'-- k-,  -  k  -ie 


Relations -  (IT)'  and  (l8  )  follow  . ..easiiy^  f^  (I?)*  :Ihus., 


lim 
G  — >  0 


.lim 
e  -» 0 

•ip*x' 


^  1  r  i£"(x“x')“  ik-x 

- q  /  - T - 5 -  Ilf (x'  )4x' <3x  dp  ' ; 

{2ny  d  k'2-  k^+ie 


®(£  “  k)t(x'  )dx'dp 


k  -  p  +  ie 


-ik-x' 


lira  /  _ 

^  J  k'2  _  k^+  ie 


i)f(x'  ,k'  )dx' . 
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Finally,  we  shall  require  the  identity 


(19);  .  J"u*(y,k)i|r_^(y,k')dy  =Ju(y,k')^tr*(y,k)cl^^ 

This  relation  follows  from  otu:  assumption  that  the  operator  V  is 
fiermitian,'  i.  e.  ; 


We  have,  in  fact, 

•  ^.u*(y,li)t(x/Is'  )<ii  ^ u*(x»i.)v(z,s)u(x,i£'  )<3id2 


'■'V  'V'  -i  )v*(x>i)u*{i,k)d^  •= j:  u(x,k'  )\ir*(2,k)ax.  ,■ 

.  .  To-  ccxnplete  the  proof  pf  (ll).we  apply  the  integrsil  equation  (14), ■■ 

using  the  right  side  in  place  of. u_^(x,lc)  and  its  complex  conjugate  in 
.place , of  .u^(x,k), to,  form'  . 

■  y'u*(a,£)u_j;(2s,is")^..=;B(^-'k’ )  ;  ' 

-ik'X  +  ik'  Ix-x'  I 

n-.'l-  ■■ 

(4rt)(2«)- 


;  .  _  1  '  /  e_ 


U-x'l- 


1  r  ii' '2  ■  i4|x-x' ) 

- - - 7”'“^  t+(x'»k)dx'clx 

(4rt)(2x)3/^  •  '  •'  + 


X  -x' 


■  1-  : 

X-x' 

-  ik 

X- 

■s'l 

(4.)=  J 

|2-S*I 

U-x' 

t*(x' >k)\lf^(y;' ,k’  )dx'dy'd^ 
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We  now  consider  the  last  term  on  the  right  of  (20)  and  apply 
(15)  to  the  functions  which  are  factors  of  and  in  its  integrand. 
We  have 


^-ik|x-x'  I  +  ik’  |x“X'  1 


U  -x'  I  |x  -x' 


ik^.(x'  ,k'  )d2'dx*di 


lira 

1 

p  giE‘(x-S')+ia*U-X 

e  ->  0 

(2«)^  •> 

>  (k^-p^-  i€)(k'^-  q^  +  : 

'  lira  . 

1 

p  |ij2*25'  -  ia*x’. 

€  -»  0 

(2rt)2'  ^ 

‘  (k^-p^-ie)(k*^-q^+i€) 

lira 

1 

e  ->  0 

(2«)3  - 

1  (k^-p^-i€)(k’^-p^+i€) 

-lira . 

1 

€  -♦  0 

(2rt)3  " 

1  (p^-k^+ie )(k'^-p^+i€ ) 

^*(x'  >k'  )d£dx’dx’ . 


If  we  apply  the  identity  (I6)  this  expression  heccmes , 


lira,  1- 
e  -»  0 


'2  ,2 


{ZiiY  ^  k  -k  +2i€  k  -p  -i€ 


.r_iE*(x'“S' )  *,  ,  V  , 

J  — 2'“2 - 


-llm 
€-»  0 


(2jt)^  '^k*^-k^+2ie  1'^ 


„ia*(y'-2')  .  ^ 

"Tg— 2 -  '|f+(x'»k)d^' . 

k  -p  +i€ 


If  we  now  use  (1?)  and  (I8)  it  hecomes 


-  l6  - 


^  °  k'^-k^+2ie 


A  “ikl^'-x'i  I 

^  J  1^'  -  x~|"  *- 


+lim  '  1 

■  ^  0  1  •  2  ,•  2  _ . 

k  -k  +2ie 


A  ik']x'-x'i  1  * 


We  now  substitute  from  (ik)  for  the  bracketed  expressions  and 


obtain  . 


■  “lim  •’  '  ^  1 

■  €  T-*  .0  ,  1 2  ,_2 


k  -k  +2ie:>'  '  (2rt) 


^  )di':  -Vj/r 


■'  '  .  /  A  ■  ik-  ‘x’  ■ 

k'^-k^+2i€  'J  .  .(2rt)  •  '  ■, 


■  '  g'^P  f  ^+C2L’>k0i|r*(2L',,lL)dx'  . 


Wow  by  (19)  the  second  and- fobrth.  terras  in  this . e;qpression  cancel 
each  other.  Hence  we  have  finally  for  the  last, terra  on  the. right  side  of 
(20)  the  expression 


"2^2" 


k  -k  +2ie  (2«) 


k*  )dy' 


V“-  /  ^:(xsk)d^-^ 


^  °  k'^-k^+2i€  ^  (2rt) 


O 
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If  we  make  use  of  (17)  and  (18)  we  see  that  the  two  terms  in  this 
expression  cancel  the  second  and  third  terms  on  the  right  of  (20), 
and  the  relation  (11 )  which  we  set  out  to  prove  now  follows. 

We  can  define  more  general  eigenfunctions  of  -A  and,  corres¬ 
pondingly,  more  general  eigenfunctions  of  -A +V  can,. be  defined.. by  .  . 
means  of  integral  eqmtions  slmilai'  to  (10.) .  .  For  exaimpley  we- can  , 
define  the  ■eigenfunction . 

■(  (21)",  /;u^^:(x>k ),.='  ■■  T  s’.'V; 

■  ■  ■■■  ■  .  of  -A.  The  function  t^q^(x>K)  will:  be  an  eigenfunction  of  ^A  if.  we 
.  assume  .that  ^^(k-' ,k).. has' the -.form  ..-’o' 

'■,;:(21a)  '  ni^(kSb)' “b+(kSli)  '''.&(k'2  -k^  - 

■  according  to  the  remarks  in'  (2;^).. 

We  can  then  define  a  corresponding  eigenfunction  u(x,k)  of  • 

.  ,  -A+V  as  the.  solution  of  the  integral  equa.ti bn  ; 

■  ;  .■(p)  '  .uCx^k)  =.lioi(x4);- 

If  we  replace  k  by  k’ in. (iO),  multiply  through,  by  m_j_(k,k' )  •. 

and  integrate  over  k',  we  see  fran  the  uniqueness  of , the  solution 
of  (22)  and  of  (10 )  that 


-  18  - 


(23)  u(x,k)  =  J"u^(x,k' )in^(k' ,k)dk' . 

Expressed  in  terms  of  abstract  operators  (23)  is 

(2^)  U  = 

wherein  U  is  the  operator  corresponding  to  the  kernel  ii(x,k)  and 
is  the  operator  corresponding  to,  the  kernel  m_^(k,k'  ),.  ' 

We  woiold  also  like  to  define  an  operator  M  .  such,  that' 

(25)  ■  u  .=  -  • 

We  can  do  this  by  showing  that  the  eigenfunction  u(2c»k)  also  satisfies 
an-  integral  equation  . 

■  , '  p  -ik|2i-x'  [ 

(26)  ■  •  .^(x,^)  ^  J  ~ 

where  UQ_(x,k)  is  an  eigenfunction  of  -A  given  liiy  . 

(27)  •  :  "Uq  (x^k)  =  •  -"•'"o/j  /■  e^-  m  (k'  ,k)dJc' . 

Here,  again,  ni_(k’,^)  must  have  the  form 

(27a)  ra_(k' ,k)  =  n_(k',k)6(k'2 

We  can  obtain  a  relation  for  m_(k',k)  by  setting  the  expression 
(26)  for  u(x,k)  equal  to  the  expression  (22).  We  have  then 
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(28)  u^Jx^k)  =  UQ^(x,k)--^  J 


ik[x-x'  -ik  x-x' 
e _ e  “  ~ 


-ik' -x 


'  fx'  x'  I  '  Is  -2'  I 

X  v(x',x"  )u(x’'#k)'3^'dx". 


We  now  multiply.  (28)  ty - TTr  and  integrate  with  respect  to  x. 

(2rt)3/2.  •  . 


According  to  (2?)  and  (23)  the  result  is 


-  ,  o, .  -ik'‘X  .  r  ik|x-x' I  '  -ikix-x'j 

C29)  B.fkss)  =m^(v,s)-Tj;  j  ^37i  ■  niTFr 


k  vCx',^”)u(x”,k)d2^'dx"da. 


It  is  clear  from  (28)  and  the  properties  of  the  Green's  functions 
'  ±ik(x-x'  )■  '  ■  .. 

^  [x  -x'  I —  "^hat  ^^^(x,]^)  is  an  eigenfunction  of  -A.  Hence>  the 

form  (27a)  of  m  (k',k)  we  have  assumed  must  be  correct. 

We  can  prove  this  more  directly  from  (29)  by  making  use  of  (17) 
and  (l8).  Ihus,  in  place  of  (29)  we  have 


-Ik'  'X 
e _ —  ^ 


(30)  m_(k',k)  =  m^(k'  ,k)  +  ^  ^  (2n)^'^^ 


-ik’.x' 
e  ~ 


/.  •  t  1  1.  lim 


(2n)' 


€  ->  0 


k“-k'>ie  k^-k’^-ie 


n  -ik' 'X' 

(k  -k  )  +e 


Y  i|f(x'  ,k)dx' 


- — F“  'l'(aSk)‘iS*^(k’2 -k^). 


(2n) 


3/2 
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where  we  have  used  the  identity 


The  form  (27a)  for  m_(k’,k}  follows  frcm  (30)  and  the  corresponding 
form  (21a)  for  m_|_(k’,k). 

If  either  one  of  the  operators  U  or  M_^  is  known  to  have  an  inverse 
•  then  (2li  ) 'implies  the  existence  of  the  inverse  of  the  other.  Then 
from  (25)  the  existence  of  the  inverse  of  follows.  If  in  ,  (24)  U  .. 

; is  set  equal  to  then 

'V;  ■  ^ 

In  this  case.  (30);  heccmes  -  ■ '  ^  • 


;  (32)  m_(k‘ ,k)  8(k'  -  k)  -  ani  J. 


.-ik*  'x'  .  ".p p  ■ .  • 


The  operator  S  defined  by  (32)  is  called  the  scattering  operator  (from., 
the  definition  employed  in  quantum  mecha^cs),  and  from  what  we /have 
just  observed  it  has  the  form 


(32a)  S  =  s(k',k)  =  cr(k' ,k)6(k'2 -  k^). 

The  forms  (21a),  (27a)  and  (32a)  imply  in  each  case,  we  recall  (cf.  ,  2.g), 
that  the  operator  commutes  with  the  Laplacian  A. 

I  ■ 

!  .  ■  ■  •  ■ 

I  ■  ’  •  ■  - 
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The  relationship  "between  the  kernel  of  S  given  in  (32)  and 
the  asymptotic  value  of  u_|_(x,k),  which  may  he  interpreted  as  the 
scattering  amplitude  for  an  incident  plane  wave,  can  he  seen  at  once 
fran  a  consideration  of  (6)  in  the  limit  of  large  x..  We  observe  that 


Ix-  x'  I  =  2x-2L*  +  x’^  ~  X 


x*x 


=  X  -  X  .X'  , 


where  x^  Is  a  imit  vector  in.  the  direction  pf  ,x.  . Then  from  (6  ), 


,3/a  >  ■  '  'V  ' 

(33):  ■  u^(i,ls)  f  y(x',s")u,(x",k)dx’dx”- 


ik-x  .  v: 

e -  (-2n)' 


'•  ikx  *  * 

, The -factor  of  ^ — -  in  the  second  term  on  the  right  of-  (33)  is  the 
scattering  amplitude  --  —  i^.(k',k),  where 


n  -ik*  ‘X*- 

(3^)  ,k)-=?  ■’^(xSx")u^(x",k)dx'dx'’ 


■and  it  is  assumed  that  k'=k  and  k'  has  the  direction  Of- x.  Then  we 
can  write  for  the  kernel  s(k‘,k)  of  S  the  relation 


(35)  -s(k'  ,k)  =  6(k' -  k) -23tit^(k'.,k)6(k'^-  k^) 


We  have  from  (25),  (32)  and  the  manner  of  deriving  the  relation 
for  S  that 


(36) 


=  U_S, 


H* 
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and  hence  that 

(37)  S  =  u\. 

We  have  used  the  fact  that  U_  is  iinitary  in  deriving  (37)>  which  in 
turn  implies  that  S  is  unitary: 

ss*  =  (u*u_^.)(u*u^)*.  =  u*(u_^u*)u_=  u*u_  =  I.  .  ■ 

We  can  conclude  from  this  that  S  is  xinitary  because  we  teve  assumed 
that  U_^  and  U'  have,  left  - and  right  inverses.  ■  • 

■  ^ora  {2k)  and  (25)  we  have 

UM  =UM.. 

Thus,  ■■ 

M_  =  ' 

from  which  we  get 

(38)  S  =  M_m“^  . 

We  can  now  state  an  integral  transform  theorem  in  terms  of  an 
arbitrary  eigenfunction  u(x,k)  of  the  operator  -A+V. 

Given  an  arbitrary  function  f(x)  of  the  type  we  consider  as 
representing  the  abstract  linear  space  (e.g.,  square  integrable  functions) 
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we  can  write 


(39)  f(x)  =f  vi(x,k)cp(k)dJlc, 


where 


(1^0) 


9(k)  =  /  fCx)  ^  /  w(k,k' )u*(x,k' )dk' >  dx 


and  the  function 


w(k,k'  )u{x,k'  )dk' 


is  an  eigenfunction  of  -A +V.  '  This  catt  be  guaranteed  by,  reqtiirihg 
that  w(k,k' )  have  the  form 


(1^1) 


)6(k2  -  k' ^ . 


It  will  also  turn  out  that  the  operator  W,  which  we  shall  refer'  to  as. 
the  spectral  weight  operator,  -corresponding  to  w(k,k^ ),  is  Herndtian  . 
and  positive  definite. 

We  can  define  W, abstractly. by 


W  ; 


then  from  (2k)  and  (25)  we  obtain 


W  = 


(43) 
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The  fact,  that  W  as  defined  by  (42)  or  (43)  is  Hermltian  and  positive 
definite  is  obvious. 

The  relations  (24)  and  (25)  also  imply 

(44)  WU*U  =  I 
and 

(45) '  uwu*  =  I.  .  . 

Relations  (44)  and  (45)  are  equivalent  to  the  transform  theorem  relations 
(39)  and  (40).  .  ' 

ii*Qm  .(42)  or  (43, )  we  see  that  the  spectral  weight  operator  W,- 
:•  needed  for  the  integral  transform  theorem'  corresponding,  to  'a  given 
eigenfunction .u(x,k)  can.be  determined, if  we  know  the  integral  equation, 
of  , the  form  (22)  which  u(x,k.)  satiafiesi  We  can  also  obtain  the  ■ 
spectral  weight,  operator  W,  in  terms  of  the .asymptotic  amplitude  of 
u(x^k)  X  becomes  large.  The  scattering  operator  S  plays  a  signi-  •  .  . 
ficant  role  in.  this  relationship,  and  we  can  therefore  obtain-  a  direct 
connection  between  S  and  W. 

Consider>  for  example,  an  eigenfunction 

(46)  .  .  ■  u(x,k)  =  y.u^(x,k' )m^(k',k)dh'  .  . 

defined  to  have  the  asymptotic  behavior 

/  ,  N  ik-x 

u(x,k)  ~  e  - 

as  X  becomes  large  in  a  "forward"  direction,  that  is,  for  the  polar  angle 
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0  of  X  in  0  <  0  <  2'  call  write  instead  of  (46),  with  ni^^(k'  ,k) 


-1 


as  the  kernel  corresponding  to  , 


(47) 


u(x,k'  )iti^^(k' ,k)dk'  =  u_^(x,k).. 


¥e  now  consider  the  asymptotic  form,  of  (47)  for  large  x,  using  (33) 
and  (34)  and  (A. l)  of  Appendix  II.  We  have,  assuming  that  in_^^(k',k) 
.has  the  form 


ii+^(k'jk)  =■  h^^(k' >4)6(k'^-  k"^). 


r2  .,2'. 


(48) 


u(x,k'  )m_^^(kf  ,k)di' 


e^~  ^(k',k)dk' 

-f*  “*  ^ 


-n 


2jt  jt  “  :  .1  t  ■ 

n  ik '  r  X 


0  0  o 


(2:t) 


m’^(k' ,k)k'^sin  0'dk'd0'dcp' 


Jti 


2rt  •  Jt  ,  .  ■  00 

dtp'  f  d9'  f  dk'  |i^.^(k' >k) 


ikx  -ikxl 

X  .(-5(k'-k)6(0'-e)6(cp*-q})  +6(k'-k)6(jt-0-0' )6(cp±jt-cp' )^ 


Jti 


(2jt) 


372 


ikx 


-n'"(.ko,k)V 


-ikx 


X 


where'  k  is  the  vector  “  x.  In  the  calculation  of  (48)  we  have  used  the 
.  -o  X  ~ 


identity 


6(k'^-k^)  =  8([k'+  k][k'  -k])  =  ~  &(k'  -k). 


0 
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We  also  have  from  (33) 


2rtl 


ikx 


(49).  u^(x,k)  ~  {-8(0-G^)5(q).(p^)  — +8(7t-e-e^)6(cp±«-9^)  ^ 


(2rt)' 


-ikx 


V-, 


./2rt>3/2  ikx 

-  tjk  ,k)  - - 

4n  o'-  X 


By  combining  (.47),  (48)  and  (4-9)  we  obtain 


n:"(k  a)  -  .  -  5(9-e  )8(cp-cp  )  -2nit  (k  >k), 

(50)  °  k  sin  G.  o  .  o  . 


/  ^:^+^(-ko■;k),  =  2^^  '•  '  G-0^)5(cp±3t-cp^) 


k  sin  0 


o 


Finally,  from  (35) »  (5p)  and  the  identity 


l  kj  =  &^k*-k)5fe'-9)5(cp’  -tp) 


k^sin  0 


we  have 


-1/ 


m_^-(k^,k)  =  5(k^- k)-2rtit_^(Js^,k)6(k^-  k^)  =  s(k^,k),  0  <  0^  <  | 


(51) 


m;^(-kQ,ls)  =  5(iSo+  ls)  or  =  &(k^3-  k),  |  <  6  <  «• 


That  is. 


=  s(k^,k)Ti  (|  -e^  +  5(ko-k)Ti  (©o-f)  • 


(52) 


o 
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Then  frorn  (42)  we  obtain  the.  kernel  corresijonding  to  W  in  the  form 
w(k^,k)  =  )m*"^(k,k'  )dk' 

=  J  s(kQ,k’  )s*(k,k‘  )<3i’4(|  e)+  s(k^,k)Ti(|  -  e2)n(e-|) 

(53) .  ■  , 

■^s"(k,k^)n(|-e)n(0^.0 

-  ^(iSo'^is)  +  s(k^,yn(|-eo)^(e-f)+s^k,k^^^^ 
wherein. we  have  used  the  unitary  property  of:  the  scattering  operator  S. 

4.  The  inverse  urohlem  • 

.In  the, inverse  problem  our  goal  is- to "obtain  the  operator  V 
Of  (2),  that  is,  the  kernel  v(x,x' ),  froiii  a  knowledge  of  the  scattering  • 
operator  S,  or  actixally  ^om  a  knowledge  of  some  part  of  S.  It  will 
be  convenient  to  rephrase  the  problem  in  a  different  form  in  terms  of 
the  spectral  weight  operator  W.  We. can  do  this  by  making  use  of  (53)>. 
and  we  observe  that  w(k^,li)  corresponding  to  W  requires  only  part  of  S 
for.  its  complete  specification. 

In  the  new  form  the  inverse  problem  is  to  find  V  from  a  knowledge 
of  W.  For  this  purpose  it  -vdll  be  convenient  to  rewrite  (44)  and  (45.) 
in  the  form 

(54)  wir^  = 

u  =  u"^  .  • 


(55) 
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If  a  solution  U  of  (54)  and  (55)  can  be  fomd  then  V  will  be 
given  by 

(56)  V  =  -U  A  +  A  . 


We  can  see  this  by  applying  -U  A  to  the  function 


ik-x' 


)  =  /  Ti(x>k)  - dk  =  U  , 


where  here  u(x,x' )  is  the  kernel  of  the  integral  operator  representation 
of  U  on  the  space,  of  functions  of  x.  By  (55)  the  result  will  be 


ik'x'  .  p., 

'  -  dx*  =  -k^ 


■ik‘x' 


.u(x,x'  )  A  - 


u(x>x'  )  — 775  =  -k^u(x,k) 

(27t)3/2 


That' is,  we  shall • have  the  expression' equivalent  to  (56): 

-Au  (x»k)+  r v(x,x*  )u(x%k')<3x'  =  k^u(x>k)  ,  ■ 

_  ■  ■  •  U  '  . 

where  ■v(x,x').is  the  kernel  corresponding  to  the  operator  Y  given  by  (56). 

From  (54)  and  (56)  it  follows,  since  A  is  Hermitian,  that  if  V 
is  Hermitian,  i.e.,'  . 


(57), 


V  =  , 


U  A  WU*  =  uw  A  u*  , 


then 
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where  we  have  used  the  fact  that  W  is  Hermitlan.  Since,  according 
■to  (55),  U  has  an  inverse  it  follows  that 

(58)  AW  =  WA  , 

which  means  that  the  kernel  corresponding  to  W  in  the  k  representation 
has  the  form 

(59)  .  w(k,k' )6(k^-  k'^) 

as  we  have  assimied. •  Relations  (58)  and  (59)  are,  in  fact,  necessary 
and  sufficient  for  (57)  to  hold. 

Obviously  the  solution  of  (5^)>  (55)  ior  U  is  not  unique  for 
.  a  given  W. '  As- we  have  Just  observed  when  W  has  the  form  (59)^  which 
we  always  assume,  then  any  V  corresponding  to  a  solution  U  will  auto¬ 
matically  be  Hermitian..'  In  order  to' guarantee  that  V  be  truly 
sejLf -ad joint  and  that  the  correct  scattering  operator  be .reproduced 
we  must  introduce  other  conditions.  The  general  question  of  the. 
uniqueness  and  existence  of  Y  luider  any  general  condition  which  leads  ■ 
to  the  correct  scattering  operator  is  still  open,  and  we  shall .not 
attempt-  to  settle  it  here.  Instead  we  shall  consider  a  particular 
condition  which  leads  to  a  tinique  'V  such  that  the  part  of  the  ' scattering 
operator  needed  to  define  W,  as  in  (53),  'vrf.ll  be  reproduced. 

The  condition  we  shall  prescribe  is  a  restriction  on  the -form  of 
the  kernel  u(x,x' ).  We  -write,  first. 


(60) 


U  =  I  +  K 
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and  then  assvime  that  K  is  represented  hy  a  kernel  of  the  form 


(6l) 


K  =  q(x>x'  h(x’  -  x) 


Then 


and 


)  =  ^(x“  x*  )  +  <l(x#x'  )^(x'  -  x) 


IK-X 

(62)  u(x,k)  =  — — rr::  + 


00  2jt  jt 


(2jt) 


3/2 


ik-x' 


X  o  o 


Q.(xj2S'  )  — TTo  x'sin  ©'dO'dcp'dx' 

(2:t)3/2  .  ■■ 


From  ( 62 )  we  •  have  f  omally 


(63) 


u(x,k)  ~ 


ik*x 


(2jr) 


for  large  x  as'  required. 


Now  from  (61)  we  see  that 


(6k) 


■  T^*  = 


K*  =  q(x»x' )  =  q.*(x’  ,x)n(x  '-x' ) 


We  shall  prove  that  the  assumption  involved  in  (60),  (61)  and 
(6k)  leads  to  the  following  theorem:* 

For  any  U,  of  the  form  given  in  (60),  (61)  and  satisfying 
the  equation  (5k), 


WU  =  U  , 
o  ’ 


The  proof  of  this  theorem  was  fvirnlshed  by  Professors  J.  B.  Keller 
and  B.  Friedman. 
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we  have  the  relation 

(65)  ^ 

l.e.f  If  U  has  an  Inverse  the  equation  (55) 

U  =  U"^ 
o 

will  be  satisfied  automatically. 

PROOF;  We  have 

(66)  UU^  =  UWU*  =  U*U*  , 
since  by  (U2)  w  is  Hermitian.  On  setting 

(67) 

where 

(68)  =  q^(x^x' )ti(x' -  x)  , 
we  have  from  (66) 

K  +  +  KK^  =  k"  +  +  KjK- 


(69) 
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In  terms  of  q  and  q  this  Is 

o 

)n(x'  -x)  +qQ(xjX'  )'n(x'-x)  + J' q(x,x"  )Ti(x"-x)q^(x",x'  )ti(x'  -x"  )dx" 

~  '3-  (x' >x)’l(x-x' )  +  q*(xix)Tl(x-x' )+  y^q*(x' ^x"  )^(x"-x'  )q*(x" ,x )ti(x-x"  )dx" 

Each  term  on  the  left  vanishes  for  x’'  <  x,  and  each  term  on  the  right 
vanishes  for  x  >  x' .  It  follows  that  both  sides  are  zero  identically. 

Thus ,  we  have  proved  that 

UU'  =  I  ■+  K  +  K  +  KK  =  I 
o  o  o 

as  desired. 

We  now  observe  that  if  we  set 
■  ■  .  W  =  +  I, 

then  eqmtion  (54)  heccmes 

OK*  +  n  +  K*  =  K 
o 

or 

(70)  Kn  +  n  +  K  =  K*  . 

In  tenns  of  the  integral  operator  kernels  we  obtain  from  (70),  if  we 
use  the  fact  that  K  and  have  the  form  given  hy  (6l), 


& 

.si-, 


'A 
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2k  It  «> 


(71) 


<l(x>2"  )‘*>(x">x'  )x"  sin  0"dx"d0"dcp"  +  (iO(x,x'  )+  q(x,x'  )  =  O^x 


For  each  fixed  x  eqiMition  (71)  is  an  integral  equation  for  q(x,x'  ) 
as  a  function  of  x'  in  the  range  x'  >  x.  If  equation  (71)  has  a  well- 
hehaved  solution  for  q(X}X' )  in  the  range  x'  >  x  we  can  form  the  kernel 


n(x,x' )  =  S(x  -  x' )  +  q(x>x' )h(x'  -  x). 

We  can  prove  that  (71)  has  a  unique  solution  hy  making  use  of 
the  positive  definite  character  of  the  spectral  weight  operator  W.  It 
follows  then  frctn  (63)  that  the  solution  we  have  given  reproduces  the 
part  of  the  scattering  operator  required  to  form  the  operator  W. 

We  can  show  that  the  solution  (assumed  to  he  bomded)  q(x^x'  ) 
of  (71)  is  unique  by  a  standard  argument  which  goes  as  follows. 

Suppose  q^(x^x'  )  is  a  second  solution  of  (71)  and  tliat  q  and  q^  are 
bo\aided.  Then  the  function 


r{x,x' )  =  q(xjx'  )  -  q^(x^x'  ) 


satisfies 


2k  K  00 


I  I  f  ©"dx’'d0"d9"  +r(x,.x'  )  =  0,x'  >  : 


'>x. 


o 


o 


X 
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Whence 

y  Ti(x"-  x)r(xjx"  )  +^(2"  ~x'  )3d.x"  =  0. 

This  has  the  operator  form 

rw  =  0. 

We  multiply  (operator  multiplication)  on  the  right  by  F*  and  apply 
to  an  arbitrary  vector  cp: 

=  0. 

The  inner  product  of  cp  with  this  is 
(cp,rWF*9)  =  (r*p,WF*cp)  ■=  0. 

Since  W  is  positive  definite  the  vanishing  of  the  quadratic  form  implies 

r*cp  =  0. 

Since  cp  is  arbitrary  then 

r*  =  0, 

and  hence 

I 

)  » 

as  required. 


Finally,  we  conclude  this  section  and  the  article  with  the  remark 
that  the  relation  (56)  for  V  reduces  to  a  simple  expression  in  the  case 
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In  which  V  is  represented  by  the  kernel  v(x,x' )•  This  expression  gives 
v(x,x' )  directly  in  terms  of  q(x,x' )  which  is  a  solution  of  the  integral 
eqmtion  (71).  JVotn  (56)  we  obtain 

V  =  (M-  KA)K^, 


or  equivalently 

■''■(xjx'  )  “  ^(q(SjS'  -  x))  -  A'  {q(x,x'  )q(x'-  x)} 

+  J  [A{q(x,x")q(x"-x)}]{q^(x",x)Ti(x-x")}dx" 

-  y'A"  [{q(x,x")q(x"-x)}]  {q^(x",x)Tl(x-  x''))dx". 

If  we  express  the  f mictions  eind  Laplacians  on  the  right  in  spherical 
coordinates  then  because  of  the  fact  that  the  derivative  of  the 
fimction  is  a  8  function  we  find  that  v(x,x'  )  is  a  sum  of  two  terms  » 
one  of  which  has  6(x'-x)  as  a  factor  and  the  other  of  which  has  tj(x'-x) 
as  a  factor  (the  terms  in  &'(x'-x)  cancel  out).  Since  V  must  be 
Hermitian  we  have 

v*(x',2)  =  v(s,x' ), 

and  thus  the  coefficient  of  the  term  in  ii(x'-x)  must  vanish.  What 
remains  is  the  following  expression  for  v(2J,2'  )i 
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(72) 


v(x,x' ) 


-2S(x  -x»  ) 

2 

X 


1. 

3k 


f  N 

2 

,  X  q(x,e,(p;  x,9',(p'  ) 


where  q(x,0,cpj  x^9'  ,cp'  )  is  obtained  from  q,(x#x'  )  by  using  spherical 
coordinates  and  setting  x'  =  x.  ¥e  recall  that  q(x^x' )  is  the  solution 
of  the  integral  equation  (71)  whose  kernel  i*j(x,x'  )  is  the  representer 
of  the  operator  a  =  w  -  I  (cf.  (53) )  . 


APPENDIX  I 


Proof  of  the  Identity 


1  e 


±ik  X  -  x' 


Ij-Jt  lx  -  x' 


lim  1 
e  0 


(2«)' 


gip* (x-x' ) 

Ti  2  . , 

k  -  p  ±ie 


dp 


We  have 


^ip. (x-x'  ) 


(2n)^^  k^-  p^±ie 


dp  = 


2n  ir  oa 

1  r  r  r  e 


(2rt)-^  w  w 
0  0  0 


ip [x-x 'I cos  ©  2 

~2 - 2 -  ®  dpdGdcp 

k  -p  ±ie 


‘ - f 


(2n)  (x-x' 


j-ip|x"x'  |_  gip|x"x' 

Ti  2~~^ 

k  -  p  ±  ie 


pdp 


p  sin  p[x-x' 


dp  = 


(2jt)  |x-x'|  o  k^-p^±ie  (2jt)  |x-x' |  k^-p^±ie 


—  f 

f-v'  I 


p  sin  p1x-x' 


dp 


CO 
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2i(2jt)  |x-x' 


ipx-x'  , 
e  '  'p  dp 


”  -Iplx  -  x' I  , 
r  e  -^1-  -  'p  dp 


o  ±le  -p)(y4^±ie  +p)  -»  (/k^±ie  i€+^ 


kn  x-x' 


i  ±ik  x-x*  1  ik  x-x* 
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4jt[x  -x* 


±ik  X  -x* 


APPENDIX  2 


Proof  of  the  asymptotic  relation 


(A.l)  e^-‘%^sin  Q  ~  -2iitk6(e -G*  )5(cp- cp*  ) 


ikx 


X  ^  00 


+2ij:k8(jt-  e  -  e*  )6(cp  ±n  -  cp’  ) 


-ikx 


We  begin  with  a  proof  of  the  relation 


±i^ 


(A.  2)  a 

'  '  a  — )  00 


lim  _±iax^  _  4 


6(x). 


Proof  of  (A.2); 
Consider 


Te-^ 

a  ->  00  I 


»,  lim  i  r  ±iax^„/  ,  lim  \  P  ±iax^„, 

f(x)dx  =  ®  f(x)dx  +  „  a  /  e  f(x)dx 


a  ^  00 


P  ±it 

4  /  ^-r-  dt  +  4 
J  ^2 
■  0 


±it 


dt 


±ir 


f(0)  =  e  ^  fid), 
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where  we  have  used  the  change  of  variable 

1  1, 

X  =  a  ^  . 

The  proof  of  (A.l)  can  be  carried  out  by  the  method  of  stationary  phase. 
We  write  first 


k*x  =  kx[cos(cp  -(p*  )sin  G  sin  0'  +  cos  0  cos  0'  ]. 

Stationary  points  where  derivatives  of  k*x  with  respect  to  0  and  cp 
vanish  are 


and 


^2 


cp  +  n, 
cp  -  Jt, 


0  <  T  <  Jt 

Jt  <  cp  <  2jt 


0. 


We  consider  separately  two  power  series  expansions  of  k-x  up  to 
quadratic  terms  about  each  of  those  stationary  points  in  the  0',9' 
plane.  The  expression  e^~'^^sin  0  can  be  regarded  as  a  sum  of  two 
terms  in  which  the  phase  k*x  of  the  exponential  is  replaced  by  the  power 
series  about  each  of  the  stationary  points  in  turn.  If  this  is  done 
and  (A. 2)  is  used  we  arrive  at 
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lim 

X  00 


e^— "■^^sin  0  = 

X  ->oo 


X 


X  k^sin  e  5(0  -©'  )6(cp  -  cp'  ) 


lim  -ikx 
e 

X  -» oo 


/_2— Jl  /.  2-  Z 

■J  X  sin  0’^  y  xk  sin  0 


sin  0'6(jt-  0  -0'  )6(cp  ±jt  -cp' 


ikx  --ikx 

=  -2iitk6(0  -0'  )6(cp  -cp'  )  - -  +  2iitk6(xr  -0-0'  )6(cp  ±jt  -  cp'  )  - 
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